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AN  APPROXIMATE  SOLUTION  TO  A  TWO-DIMENSIONAL 
STATIONARY  HEAT-CONDUCTION  PROBLEM  AND  ITS 
USE  IN  THE  DETERMINATION  OP  THE  EFFECT  OF 
ARTIFICIAL  COOLINQ  OP  TURBINE  COMPONENTS 

D.  A.  Pereverzev,  Eng. 

Hydraulic  Machine  Laboratory,  Acad.  Scl.  Ukr.  SSR. 


Certain  getierallzed  solutions  to  the  problen 
of  a  stationary  tenqperature  field  In  a  semlbounded 
body  with  Internal  cylindrical  heat  sources  and 
approximate  solutions  to  problems  with  boxuidary 
conditions  of  the  third  type  are  given.  The 
equations  derived  were  used  to  calculate  the 
effect  of  cooling  a  turbine  dnan  rotor. 


Let  us  consider  the  problem  of  a  stationary  temperature  field  in 
a  semibounded  body  with  Internal  cylindrical  heat  sources. 

Let  there  be  distributed  within  the  semlbounded  body  at  one  and 
the  same  depth  h,  but  with  different  distances  between  their  axes 
(2x^  and  S;  x^  <  an  infinite  series  of  cylindrical  channels  of 

diameter  d  >  2r  (Fig.  i).  The  surface  teiqperatures  of  the  body  tp  and 
the  cylindrical  channels  t^  (boundary  conditions  of  the  first  type) 
are  given,  with  tp  >  t^  (cooling  of  the  body  by  a  system  of  submerged 
channels.  The  heat  flux  £  Issuing  from  a  \mit  length  of  the  channel 
and  the  heat-conduction  coefficient  of  the  body  X  are  assumed  to  be 


PTD-TT-62-1758/l+a44 


1' 


constant.  We  shall  assiame  that  the  axes  of  the  channels  are  parallel 
to  the  surface  of  the  body.  With  these  stipulations,  the  process  of 
heat  propagation  for  a  stationary  state  Is  described  by  the  differential 
equation  of  heat  conduction: 


The  x-2Lxls  Is  located  on  the  surface  of  the  body,  idille  the  y-axls 
Is  symmetrical  with  respect  to  adjacent  channels  (Fig.  1).  It  Is 
necessary  to  find  a  solution  to  Eq.  (l)  which  satisfies  the  given 
boundary  conditions  (tp  and  t^). 


Fig.  1.  Schematic  for  the  derivation  of  the 
temperature-field  equation. 


We  shall  use  the  method  of  superposition.  In  order  to  solve 
£  q.  (l)  .  We  shall  assume  that  the  thermal  capacity  of  the  first 
channel  Is  concentrated  along  a  line  having  the  coordinates  x  <■  x^, 
y  »  h^  and  parallel  to  the  axis  of  the  channel.  Sy  analogy,  taking 
the  condition  of  synnetry  Into  account,  we  do  the  same  thing  with  the 
remaining  channels.  We  obtain  a  system  of  point  sinks  (-q).  Synnetrl- 
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cal  with  them  we  set  up  a  system  of  sources  of  the  same  capacity  (-fq). 
The  temperature  ^  at  an  arbitrary  point  M(x,  y)  is  determined  from  the 
superposition  of  the  temperature  fields  produced  by  the  systan  of  point 
sources  and  sinks. 

We  Introduce  the  notation  9  «  tp  -  t.  Then  the  temperature  field 
of  a  sink  Is  determined  from  the  expression: 


and  the  temperature  field  of  a  sovirce  by  the  dependence 


The  total  effect  of  the  source  and  the  sink  gives  the  temperature 
field: 

1.^. 


As  a  result  of  the  combined  effect  of  all  sources  and  sinks  (Pig.  1). 
we  obtain  the  temperature  field: 


Determining  the  values  of  the  radii  from  geometrical  considerations 
(Fig.  l)  and  substituting  In  (2),  we  obtain,  after  transformations: 


/r 


In 


IS  <. + {X-  xJT  +  (*.  “if + 


+|jln 


!£jlrJ£rii!£±<5i±i?l  j. 


t5  .t+(jr  +  jr.)|»  +  (*.  +  rt» . 


+2  In 


IS 


«  I 


(?) 


PTD-TT-62-1758/1+2+4 


Here  the  first  bracket  Includes  channels  with  subscripts  ^  while 
the  second  bracket  Includes  channels  with  subscripts  11  to  the  right 
and  left  of  the  y-axls,  respectively,  (Plg.  1). 

According  to  page  282  of  Ref.  [1],  we  have 

^  *  "f"  (2te  — I  * 

Ri^resslon  (3)  Is  easily  reduced  to  this  form  by  the  substitution: 

+  = 

We  shall  convert  the  series  in  the  brackets  of  equality  (3)  Into 
infinite  products.  After  substitution  and  transformations  we  obtain: 

for  the  first  bracket 


r.  •?  If  1 

II  id  3  r 


•4  ^  4  T. - 

I +-pr  — •  ['  +  ■0inr+ v,jk  J['  +(isrrpy 


-In  cfcw.-cwV,  ^ 

Sk  ^ 

e* -J- (*»  +  *)— eo* -3- (4r  —  jr,J 

2ii^  '  ~  " 

ir)  — cot  -j-ijr  — 


analogously  for  the  second  bracket 

^2*  2b 

efc  "y  (*•+*)— ccj (jr  +  jrJ 

Si  ' -  '  '  •  t-m. 

The  heat  flux  Issuing  frcmi  a  unit  length  of  the  channel  is: 


«  ch  >  cosh 


where  Is  the  dimensionless  thermal  resistance  of  the  channel. 

0 

Substituting  (?a),  (3b )«  and  (3c)  Into  (3)  we  obtain*  after  simple 
transformations 


ch-j-  — — Xj 

ch  — cot-j-(x— xj 


X 


^  7m 

ch  -j-C*.  +  v)  —  CM  7  (x  +  X.) 

X  — 5;j - -  - 

ch  -=-  (*,  —  u)  —  cot  -s-  (X  +  xJ 
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(M 


where  T  is  the  relative  excess  temperature  at  the  point  of  the  body 
under  consideration. 

The  value  of  ^  determined  from  (4)  satisfies  eqxiatlon  (4)  and  a 
boundary  condition  of  the  first  type  on  the  surface  of  the  body  (for 
y  *  0*  T  «  1*  or  t  ■  tp). 

It  can  be  shown  that  for  »  0  and  x^  *  the  equality  (4) 
leads  to  the  solutions  obtained  In  Refs.  [2  Euid  3l  for  equally  distri¬ 
buted  channels. 

In  order  to  satisfy  the  boundary  condition  T  ■  0  (t  *  t^)  over  the 
entire  contour  of  the  channel,  we  have  at  our  disposal  In  Bq. 

(4)  only  three  constants;  x^,  h^,  and  R^.  With  their  aid  the  condition 
T  »  0  may  be  satisfied  at  only  three  arbitrary  points  of  the  original 
contour.  The  remaining  points  of  the  Isotherm  T  «  0  are  determined 
from  the  equation 


cIi-^(*,^d)  — cot-j-  (x  — X.) 

~5b - sr - 

X 

cfc^(*,  +  jr)-co«^(x  +  xJ 

^ — K - Si - 

**  X  ~  X 
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which/  In  the  general  case/  is  not  an  equation  of  a  circle. 

Therefore  Eq.  (4)  solves  the  given  problem  {approximately. 

The  approximateness  Is  due  to  the  fact  that  the  thermal  capacity  of 
a  circular  cooling  contour  Is ^concentrated  In  one  discharge.  The  saise 
thing  applies  to  the  solutions  obtained  In  Refs.  [2  and  3]. 

Thus,  Eq.  (4)  Is  an  exact  solution  of  the  problem  for  the 
contour  described  by  Eq.  (3)*  which,  however,  given  a  proper  choice 
of  the  determining  points  and  given  small  ratios  r/S,  may  turn  out  to 
be  very  close  to  a  circle.  Calculations  Indicate  that  If  the  abscissa 
of  a  sink  Is  given,  assuming  that  »  x^,  while  points  1  and  2,  located 
on  the  vertical  diameter  of  the  original  contour  (Fig.  1)  are  taken 
as  the  determining  points,  then  such  a  choice  ensures,  for  r/s  ^  0.06- 
-0.07  over  the  entire  possible  range,  a  variation  of  h/s  (h  >  r)  and 

^  i  3/4)  and  an  almost  complete  coincidence  between  a  circu¬ 
lar  contour  and  the  contour  described  by  Eq.  (3). 

Taking  Into  account  what  has  been  said,  let  us  detemJjie  the  sink 
parameters  In  Eq.  (4).  Substituting  the  coordinates  of  the  deter¬ 
mining  points  1  and  2  (Fig.  l)  and  taking  Into  accoxmt  the  properties 
of  hyperbolic  functions,  we  obtain  the  following  system  of  two  equations 

with  two  unknowns  (h^  and  R^); 

o  o 

p*int  /  x=x„  yssh  —  r,  Tt=0 


if  **<^(*.+ 

In  - - 

I .  s  (^.  ~ 


h-r) 


*  +  0 


'  — r)— 

-  mi™-  1-1  I 

ch -y  (A.  —  A  +  r) — CM ‘y  X| 


I- 


(6) 


point  2  XasJC,,  IfssA-fr.  TmtO 
•n -J- (*.  +  *  + d  Teh -J  (*»  +  *  + d  "-eo*  y  X,' 
*hy(*  +  e— chy(*  +  r-*J-CMyx, 
-2«|. 

•'■h*  Hah 
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(7) 


An  annlytlcal  solution  of  the  systea  of  equations  (6),  (7)  In 
general  form  Is  liqposslble.  But  with  •  0  and  »  S/^«  the  system 

Is  simplified.  Equating  (6)  and  (?)  and  applying  the  formulas  of 
hyperbolic  trigonometry^  after  transformations,  we  havet 


1.  With  x^ 


or 


Here  the  positive  root  gives  the  ordinate  of  the  sink,  idille  the 
negative  root  give  the  ordinate  of  the  source.  We  are  Interested 
only  In  the  values  y  >  0.  Therefoire,  we  assume  finally  that: 


Substituting  (8)  In  (6)  or  (7)»  we  obtain: 


2.  With 


S/^ 


.  *T'  * 

ll 

[efc-Jf  r  1 

^  *r’~ 

/■-T* 

“li - ' 

(8) 


(9) 


(10) 


(11) 
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For  0  <  <  s/^«  lies  between  the  values  determined  according 

to  (8)  and  (10).  This  makes  possible  a  relatively  simple  solution  of 
Eqs.  (6)  and  (7)  In  nxunerlcal  form  by  selection  of  the  quantity  h^. 
For  practical  calculations,  h^  may  be  determined  adequately  enough 
from  (10).  Moreover,  the  quantities  R^,  calculated  according  to  (6) 
and  (7),  turn  out  to  be  sufficiently  close  to  each  other.  If  they  are 
different,  then  to  increase  the  accuracy  of  the  subsequent  calculations 
of  tenpe rat\ire  according  to  £q.  (4),  it  is  necessary  to  Introduce 

their  average  value. 

If  r  is  substantially  less  than  h  and  S,  then  It  Is  possible  to 
set  h^  »  h,  slnh  ir/s  r  «  r,  cosh  r  «  1,  while  from  expression  (6) 

O  Dm 

or  (7)  we  obtain: 


|— i — I  — ; - ii - 1- 

i\  -j-f  J  l-cot-y-or,  J 


(12) 


The  degree  of  accxiracy  of  the  solution  may  be  Increased,  If  It 
Is  constructed  with  the  aid  of  several  heat  sinks.  Let  the  number  of 
sinks  be  equal  to  m.  In  this  case  there  will  be  ^  constants  which 
make  It  possible  to  satisfy  the  boundary  condition  T  «  0  at  ?m  points 
of  the  given  circle.  It  Is  possible  by  such  a  method  to  obtain  a 
solution  of  the  problem  for  channels  with  contours  of  arbitrary  shape. 

Applying  the  method  of  superposition  for  an  arbitrary  point  N 
(Fig.  2),  we  obtain  in  this  case: 


"  cb -7  (**<+»)- cos -T-U-jr.!) 

g - 


(15) 


or  In  relative  units 


-Sacte— JJ-- - B - . 


(1^) 


where  m,  and  are,  respectively,  the  number,  coordinates, 

and  dimensionless  thermal  resistances  of  the  sinks  in  the  contour;  S 
is  the  spacing  of  sinks  of  equal  depth  along  the  x-axls; 
is  the  productivity  of  the  sinks. 


The  total  thermal  resistance 
of  a  channel  with  a  contour  of 


arbitrary  shape  ist 
I 


(15) 


Pig.  2.  Cooling  of  a  semibounded 
by  a  system  of  channels  with 


Here  the  y-axls  is  directed 


contoiirs  of  arbitrary  shape. 


through  the  center  of  gravity  of 
the  transverse  cross  section  of 

the  chaimel,  but  it  may  be  directed  arbitrarily.  The  argumentation  la 
not  changed  by  this. 

In  a  more  general  case  there  may  be  several  series  of  channels 
which  differ  in  the  shape  of  their  contour  and  in  their  spacing.  In 
this  case  2uldltlonal  sums  will  appear  in  £q.  (15)« 

Every  sink  is  characterized  in  (i4)  by  three  parameters: 
h^^  and  R^^*  which  are  still  unknown.  In  order  to  determine  them, 
it  is  necessary  to  set  up  a  system  of  equations  from  the  condition 
t  *  0  with  respect  to  the  individual  points  of  the  cooling  contows. 

Given  the  number  of  sinks,  it  is  possible  to  include  any  desired 
number  of  points  of  the  contoxu*.  Put  in  this  way,  the  problem  may  be 
solved  with  any  desired  degree  of  accuracy  necessary  in  practice.  Itae 
difficulties  of  solution  of  complex  systems  of  algebraic  equations  in 
the  determination  of  ^m  unknowns  may  be  eliminated  by  the  use  of 
electronic  computers.  In  solving  the  system,  it  is  necessary  to  take 
into  account  that  the  heat  sinks  must  be  distributed  within  the  r^ion 
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enconpasaed  by  the  original  contour. 

In  a  number  of  cases  the  solution  of  the  system  Is  simplified. 

Thus,  for  example,  for  symmetrical  contours,  the  number  of  unknoNns 

may  be  reduced  by  half.  The  number  of  unknowns  may  also  be  reduced. 

If  the  abscissas  of  several  sinks  are  given.  For  circular  contours. 

if  r^  Is  significantly  less  than  h^^  and  It  Is  possible  to  set 

h^^  -  h^  (the  sink  coincides  with  the  center  of  the  circle),  and  then 

everything  reduces  to  a  solution  of  a  system  of  linear  algebraic 

equations  with  the  xmknowns  z-  ■  . 

*  Rol 

Let  us  note  some  of  the  most  Important  properties  of  the  tempera' 
ture  field  described  by  Eq.  (1^). 

1.  The  tenperature  of  a  body  at  Infinity  Is  determined  by  the 
relationship 


im  ^ 

—  —  .Jl 


(16) 


2.  Let  us  consider,  for  any  below  the  cooling  channels,  the 
heat  flux  and  the  average  temperature  of  the  body  In  a  step-type  cross 
section.  In  the  general  case,  their  values  are  expressed  by  the 
relationships! 


heat  flux 


0 - ^(fr  — 

(17) 

average  temperature  of  the  body 

t 

7.=-ffrdx. 

(18) 

f 

Then  , 

1  f 

(18«) 

But  below  the  cooling  channels  for  any  Q  ■  0.  and  therefore 
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^1^  ■  0,  or  *  const  •*  «  .  ¥lth  a  reduction  In  ^  from  the 

value  (y^^  Is  the  ordinate  of  the  lower  point  of  the  deepest  contoitr), 
r.  starts  to  increase.  Consequently,  below  the  cooling  channels  the 
average  temperature  of  the  object  In  a  step-type  cross  section  (desig¬ 
nated  by  has  a  mlnlmun  value  and  Is  equal  In  all  cross  sections 

to  the  temperature  at  Infinity. 

3.*  The  nonunlfomlty  of  the  temperature  with  respect  to  the 
spacing  decreases  rapidly  with  increasing  distance  from  the  cooling 
channel  into  the  interior  of  the  body.  Calculations  Indicate  that  at 
the  relative  depth  J  i  +  (0*^  “  0-5)  (h^j^  Is  the  ordinate  of  the 
lower  sink)  the  temperature  field  at  all  points  of  the  body  Is  practi¬ 
cally  leveled  to  the  te^perat\ire  at  Infinity. 

JBy  virtue  of  the  above-enumerated  properties,  the  temperature 
ty^  »  Tg  may  serve  as  a  practical  criterion  of  the  effectiveness 
of  cooling  a  semlbounded  body  with  a  system  of  submerged  chaxinels.  The 
relatively  simple  expression  (l6)  significantly  simplifies  calculations 
In  the  choice  of  the  optimum  dimensions  of  the  channels.  This  tempera¬ 
ture  may  serve  as  a  criterion  for  the  most  rational  choice  of  the 
number  of  sinks.  For  example.  If  the  furthest  Increase  In  the  number 
of  sinks  does  not  alter  essentially  the  quantity  t^  then  It  may 
be  asstuned  that  a  satisfactorily  accurate  solution  has  been  attained. 

A  third  propa>ty  of  the  field  makes  It  possible  to  solve  (l4)  for 
the  calculation  of  teB|>erature8  In  flat  walls  with  a  thickness 
^  ^ +  (0. 4-0.5),  If  the  lower  surface  of  the  wall  opposite  the 

heat-supplying  side  Is  free  from  heat  exchange. 

If  the  temperatures  of  both  surfaces  of  the  wall  tp^  and  tpg  a^re 
given,  where  t^2  ^  either  larger  or  ssutller  than  t^,  then  with 
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^  ^ +  (0,4-.0.5)  the  equation  of  the  tetqperature  field  nay  be 

expressed  accurately  enovigh  In  the  form  of  a  sum  of  partlculau?  solu¬ 
tion:  /ri— fa 


ch  -J-  (*„  4  yj  _  cot  -|r  (jr  —  jr,i) 

55  S - 

ch  -j-  (*„  —  jf)  -  cot  Xtt) 


(19) 


The  solution  of  (19)  satisfies  Eq.  (1)  and  the  boundary 
condition  on  the  upper  surface  of  the  wall,  when  y  «  0  and  t  >  tp^. 
To  determine  the  sink  parameters  and  the  constant  C,  It  la  necessary 
to  use  the  boundary  conditions  t  »  t^  along  the  cooling  contour  with 
y  s  H  and  t  -  tpg*  latter  condition,  taking  Into  account  the 

first  and  third  properties  of  the  temperature  field  described  by 
Eqs.  (13)  and  (1^),  may  be  written  thus: 
for  y  «  H 


^.-fr.=  2]^+C//. 

fm»| 


(20) 


For  a  semlbounded  body  and  flat  wall,  a  comparison  was  made  of  the 
circle  and  the  curved  Isotherms  t  =  t^  obtained  from  Eqs;  (13) 
and  (19)  for  m  »  1  and  x^  >  0  (equally  submerged  and  equally  distributed 
cooling  channels).  As  can  be  seen  from  Pig.  4,  the  Isotherm  t  »  t^, 
with  a  reduction  In  r/s.  In  both  variations  approximates  a  circle  more 
eind  more.  With  r/s  «  0.03,  all  three  curves  practically  merge  Into 
one.  On  the  whole,  for  r/s  ^0.1  the  replacement  of  the  thermal  capac¬ 
ity  of  a  circular  channel  with  a  single  point  heat  sink  ensxires  good  accuracy 
of  solution.  Alternative  calculations  Indicate  that  this  holds  true 
even  for  h/S  <0.3,  provided  that  the  condition  h  >  r  is  observed. 

For  h/s  >  0.3  a  good  accuracy  of  solution  may  be  obtained  even  for 
r/s  >  0.1. 
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Fig.  Degree  of  approximation  of  the  Isotherm 
t  -  t^  to  a  circular  contour  (equally  submerged 

and  equally  distributed  channels)  In  the  presence 
of  a  single  heat  sink  In  the  contour. 

I)  In  a  semlbounded  body;  II)  circle;  III)  In 
a  flat  wall  (H  »  2h,  tp^^  “  ^2  “ 

1,  2)  analogous  to  Fig.  1;  5,  4)  location  of 
point  heat  sinks  In  a  semlbounded  bo<^  and  a  flat 
wall  respectively.  Initial  data:  h/s  ■  0.5j, 

WS  »  1;  a  -  r/s  -  0,2;  b  -  r/s  »  0.1;  c  -  r/S  » 

-  0.05. 


In  practlce^we  encounter  more  frequently  boundary  conditions  of 
the  third  type,  where  the  temperatures  of  the  media  adjacent  to  the 
heat -supplying  (tp)  and  heat-removing  (t^)  surfaces  are  given. 

Below  Is  a  proposed  solution  of  such  problems  by  the  method  of 
the  additional  wall.  It  consists  of  the  following. 

A  layer  of  metal  Is  added  arbitrarily  to  the  heat-supplying  side 

The  origin  of  the  coordinates  Is  raised  by  the  amount  A,  and  the 
depth  of  the  channel  must  now  be  taken  as  hp  >  h  -f  A.  For  a  flat  wall 
In  the  presence  of  boundary  conditions  of  the  third  type,  layers  of 
metal  are  added  to  both  surfaces  of  the  wall: 
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An  additional  cylindrical  layer  Is  arbitrarily  added  to  the  heat- 
removing  side,  ^e  Inner  diameter  of  this  layer  la  determined  by  the 
relationship 


3cX(/.  —  j  ^  ^ 

d — — ('■  W. 


whence  we  obtain: 


(23) 


In  formulas  (2i)  and  (^) 


(the  Blot  criteria  for  the  heat-supplying  and  heat-removing  sides, 
respectively); 

X  «  the  heat-conduction  coefficient  of  the  body  determined 
from  the  average  tenperature  of  the  body; 

Op  and  a  =  the  heat-transfer  coefficients  on  the  heat -supplying  and 
^  °  heat-removing  sides  (taken  as  constants); 

tj.  «  the  average  tanperatxure  of  the  surface  of  the  cooling 
“  channel; 

S  the  spacing  of  the  cheumels; 

d  >  the  diameter  of  the  cooling  channel; 

dp  =  the  hypothetical  diameter  of  the  cooling  channel  (innor 
^  diameter  of  the  additional  cylindrical  layer) 

A  channel  with  a  contour  of  arbitrary  shape  may  be  replaced  by 

a  system  of  circular  channels  inscribed  In  It,  provided  that  the  total 

perlaeter  of  the  circular  channels  and  the  heat-transfer  coefficient  in 

them  axe  kept  equal  to  the  original  (u  and  a^).  In  this  case  the 

hypothetical  diameter  Is  determined  by  the  formula 

where  d.  -  diameter  of  one  of  the  circles  replacing  the  original 
“  contour 


14 


<■  Biot  criterion  on  the  heat-removing  aide 

®  nr" 

Jlv 

d-  ■  ~  ■  effective  dJa  meter  of  the  original  contour 
o  u 

>  heat-transfer  coefficient  in  the  original  channel 
In  formulas  (4)  and  (20),  h  and  d  are  replaced  by  h^  and  d^,  and 
the  problem  is  thus  reduced  to  a  problem  with  boundary  conditions  of  , 
the  first  type.  Consequently,  for  boundary  conditions  of  the  third 
type,  T  no  longer  depends  only  on  geometrical  factory  but  also  on  the 
magnitudes  of  Bip  and  Bi^. 

Analysis  of  alternative  calculations  shows  that  with  the  substi¬ 
tutions  Indicated,  a  small  variation  in  the  thermal  resistance  of  the 
body  Itself  occurs.  This  is  expressed  in  a  certain  deformation  (as 
compared  with  the  real  pictiare)  of  the  temperature  field  near  the 
original  channel,  which,  for  all  practical  piirposes,  has  little  effect 
on  the  accuracy  of  the  calculations.  With  an  Increase  in  the  external 
resistances  (with  a  reduction  of  Op  and  especially  a^)  the  influence 
of  these  errors  decreases,  and  the  accuracy  of  the  calculations  increases. 
The  same  thing  takes  place  with  an  increase  in  the  depth  of  the  channels. 
Thxis,  even  with  h/s  >  0.|^  a  contour  of  arbitrary  shape  may  be  replaced 
by  a  single  circular  one  at  the  same  depth  h  without  great  loss  of 
accuracy  in  the  solution  (the  determination  of  the  cooling  effect). 

Since  there  is  no  strict  solution  to  the  problem,  the  limits  of 
applicability  of  the  method  of  the  additional  wall  can  be  established 
only  through  comparison  with  experiments  and  with  other  more  accurate 
methods  of  calculation  (e.g.,  the  ETA  method). 

The  solutions  derived  were  used  in  the  calculation  of  the  cooling 
effect  in  a  drum  rotor  of  a  steam  turbine.  — 

Figure  5  shows  an  element  of  the  cooled  drum  rotor  of  the  SKR-iOO 
XhTdZ  turbine  designed  for  the  initial  steam  parameters  of  y>0  atm  (abs) 
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and  630°  C.  Cooling  la  aecoqpllahed  by  a  stream  of  steam  of  reduced 
temperature  (525^  C)  flowing  along  specially  provided  channels  In  the 
tsdls  of  working  blades  1  and  In  Intermediate  Inserts  2,  which  shield 
the  surface  of  the  rotor  In  the  zone  of  distribution  of  guide  vanes  3 
[  3  J*  A  section  of  an  element  indicated  by  dot-daSh  lines  and 

solid  shading  was  studied  by  the  ETA  method  on  solid  and  plane  models 
In  electrolytic  baths  [3a* 

For  the  solid  model  It  was  assumed  that  the  drum  rotor  was  made  of 
perlltlc  steel,  and  the  blades  and  Inserts  of  austenitic  steel.  For 
the  plane  model  (semlbounded  body)  the  heat-conduction  coefficient  X 
was  taken  equal  to  the  heat-conduction  coefficient  of  the  material  of 
the  austenitic  blades.  The  depth,  the  dimensions  and  also  of  the 
silt  In  the  plane  model  were  taken  the  same  as  for  the  working  blades 
In  the  solid  model.  The  spacing  was  taken  equal  to  the  root  spacing 
of  the  working  blades. 

Figures  3a  and  3h  give  the  temperature  fields  obtained  on  solid 
and  plane  models  for  nominal  working  conditions  on  a  section  of  the 
second  reactive  stage.  Comparison  of  them  shows  that  below  the  cooling 
channels  the  tei!Q>erature  of  both  models  levels  off  rapidly,  while  the 
average  temperature  Biln^  plane  model  Is  only  0.01-0.2 

higher  than  the  temperatinre  of  the  rotor  (O.I73  as  against  0.13-0.16 
This  makes  It  possible.  In  the  given  case,  to  substitute  a  plane 
problem  for  a  solid  one  and  to  consider  the  temperature  of  the  rotor 
as  the  temperature  of  a  semlbounded  body  at  Infinity. 

«  It  was  asstsned  that  on  the  boundaries  of  an  Isolated  elejamt 
axial  temperature  gradients  do  not  exist.  Experiments  [3]  have  shown 
that  the  temperature  of  the  rotor  drum  varies  Insignificantly  from 
stage  to  stage  (it  decreases  with  the  flow  rate  of  the  steam),  so 
that  this  assumption  was  Justified. 
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Fig.  Element  of  a  cooled  rotor  studied  by  the 
El  A  method. 

a)  temperature  field  In  the  cross  section  EP 
obtained  on  the  solid  model;  cooling  channels; 

In  the  blade)  plane  silt  6  x  16  mm;  In  the  lixsert) 
circular  channel  with  3  ram  diameter;  b) 
temperature  field  obtained  on  the  plane  model, 
B1j,/b1q  «  19.15/2.25;  cooling  channel)  plane  silt 

6  X  16  mm. 


In  accordance  with  (6),  (10),  and  (I6),  relative  temperatures  were 
calculated  for  the  nominal  and  certain  variable  regimes.  The  plane 
silt  was  replaced  by  two  clrculauc*  chemnels  while  keeping  the  original 
u  and  a^.  The  basic  geometrical  dimensions  of  the  slit  and  the  results 
of  the  calculation  eu*e  presented  In  the  table.  Also  presented  In  the 
table  are  ETA  data,  which  were  obtained  for  analogous  Blp  and  Bl^  at 
the  Department  of  Turbine  Construction  of  the  Kharkov  Folytechnleal 
Institute  by  V.  M.  Kaplnos  and  Ye.  I.  Bublikov.  The  results  of  the 
calculations  and  the  ETA  data  coincide  almost  completely. 

It  should  be  emphasized  that  the  plane  model  under  consideration 
takes  Into  account  well  enough  coqplex  heat  exchange  at  the  root  of 
the  blades  and  In  the  transition  chamber  4  (Fig.  5)  for  the  acttial 
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case  tdilch  was  Investigated  In  the  literature  [3].  However^  it  is  not 
Ijqpossible  that  for  other  geometries  of  the  cooling  channels  and  other 


boundary  conditions,  the  heat-transfer  results  obtained  on  plane  and 


solid  models  will  differ  more  significantly. 
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Nevertheless,  the  computational  method  proposed  may  be  recoomended 
for  the  selection  of  the  optimum  dimensions  of  cooling  channels  and 
for  estimations  of  the  temperature  of  a  cooled  drum  rotor  In  designing 
new  hlgh-temperature  turbines,  if  a  system  of  cooling  amlogous  to  the 
one  under  consideration  is  used. 

The  solutions  obtained  may  also  be  used  for  calculating  the  tempera¬ 
ture  of  disks  cooled  by  passing  an  appropriate  medium  through  the  flow 
areas  of  the  tall  sections  of  the  blades  or  Immediately  on  the  rim  of 
the  disk  (assembly  clearances,  borings,  etc.).  Moreover,  if  there  is 
an  absence  of  additional  heat  removal  from  the  web  of  the  disk,  then 
solutions  of  type  (4)  and  (l4)  are  applicable.  In  the  presence  of 
additional  heat  removal  from  the  web  of  the  disk,  the  equations  of  the 
tenperature  field  may  be  obteiined  by  "stitching  together”  solution  (19) 
and  the  widely  known  solution  of  the  one-dimensional  problem  for  a 
disk  [6].  The  first  of  them  will  give  the  temperature  distributicm 

in  the  zone  of  the  rim,  the  second  will  give  the  temperature  distri¬ 
bution  over  the  web  of  the  disk. 
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Finally,  Eqi.  (1^)  and  (19)  permit  determination  of  the 
thermal  resistances  of  the  channels  (or  tubes)  with  a  coiqplex  cross- 
sectional  profile.  If  these  channels  are  laid  in  a  semibounded  body 
or  a  plane  wall. 
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OPTIMUM  DIFFUSER  SHAPE 


I.  Ye.  Vlrozub  and  A.  Sh.  Dorfman 

It  Is  well  known  that  the  basic  type  of  loss  In  diffusers  Is 
expansion  losses.  For  the  case  where  the  diffuser  has  stralght-llne 
generatrices  (flat  and  conical)  and  the  fluid  is  incompressible, 
expansion  losses  are  usually  determined  by  shock  softening  factor 
9  The  essence  of  this  method  of  determining  losses  Is  that 

expansion  losses  are  equated  with  shock  losses  determined  by 

the  Borda-Carnot  theorem: 

(1) 

where  p  Is  the  density;  and  stream  velocities  on  Intake  and 
exit. 

It  Is  considered  that  expansion  losses  In  diffusers  with  recti¬ 
linear  generatrices  make  up  part  of  the  shock  losses: 

(2) 

Many  experiments  (at  least  for  an  Incompressible  fluid)  have 
ascertained  that  factor  9  in  Eq.  2  chiefly  depends  on  the  diffuser's 
expansion  angle  6L  Figure  1  shows  the  relationship  f  (9) . 

No  doubts  arise  that  the  magnitude  of  the  losses  in  the  diffuser 
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In  Its  prescribed  length  and  degree  of  expansion  will  depend  only  on 
the  shape  of  the  walls  forming  this  diffuser.  Therefore  the  question 
comes  up  as  to  the  selection  of  the  medium  for  all  diffusers  of  a 
prescribed  length  and  degree  of  expansion,  the  expansion  being  such 
that  losses  In  It  are  the  least.  This  problem  has  attracted  the 
attention  of  many  researches.  In  Idel'chik  [2],  for  example,  thes 
is  a  survey  of  research  on  determining  diffuser  losses;  and  diffusers 
with  stralght-llne  and  curvilinear  generatrices  are  compared,  as  welT 
as  diffusers  with  abrupt  changes  in  cross-sectional  area.  But  the 
problem  of  choosing  the  diffuser  shape  In  which  losses  are  minimum 
Is  not  posed  In  this  research. 

Formula  (2)  must  not  be  directly  used  In  computing  a  diffuser 
with  curvilinear  generatrices.  Therefore  In  the  following,  when 
solving  the  problem  of  optimum  diffuser  shape,  we  use  the  concept  of 
local  expansion  angle  Introduced  In  Borfman  et  al.  (4]t 


where  F  Is  the  area  of  the  longitudinal  cross-section;  U,  the  wet 
perimeter;  and  x,  the  longitudinal  coordinate. 

It  Is  easy  to  show  that  for  an  axlsymmetrlcal  diffuser  the 
foregoing  expression  leads  to  the  forms 

_  •  (  3) 


where  R  Is  the  radius  of  the  cross-section  of  the  diffuser. 

We  will  examine  an  axlsymsmtrlcal  diffuser  of  prescribed  length 
and  degree  of  expansion  ( Fig.  2) .  The  shape  of  Its  generatrices 
must  be  detenslned  from  the  condition  that  expansion  losses  have  a 
minimum.  The  fundaswntal  region  of  a  diffuser  of  length  dx  may  be 
seen  Isolated  In  Fig.  2  regarded  as  a  conical  diffuser  with  Its 
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Pig.  1. 


Pig.  2. 


g.  5 


expansion  angle  determined  by  Formula  (?).  Differentiating  Eq.  (2) 
taking  ( l)  Into  consideration  we  will  find  that  expansion  losses  In 
this  region  amovuit  to  the  following  (more  detail  In  Dorfman  et  al. 

m): 


where  factor  9  Is  determined  by  the  value  of  the  local  expansion  angle. 
The  stream  discontinuity  equation  for  an  axl symmetrical  diffuser  may 
be  presented  In  the  formt 


(5) 


We  refer  all  the  variables  having  length  to  the  radius  of  the  entry 
section  of  the  diffuser* 

In  the  following  we  will  drop  the  strokes  above  the  dimensionless 
magnitudes. 

From  ( 4)  and  ( 5)  It  follows  that 


Therefore  loss  factor  c— f- ^  will  be  computed  according  to  the 
formula 


t-4 


(6) 


Thus,  the  problem  of  finding  the  optimum  diffuser  shape  reduces 
to  finding  a  cux*ve  R  >  r(]0  such  that,  when  Integrating  along  It, 
Integral  (6)  will  assvune  Its  minimum  value,  l.e.,  towards  the  varla> 
tlonal  problem. 

It  Is  well  known  (e.g. ,  cp.  Lavrent'ev  and  lyustemlk  [5])  that 


the  Integral 


assumes  Its  extreise  value  on  the  y  -  y(x)  curve  on  condition  that 

(d/dx)f^  -  -  0  (Euler's  equation)* 

where  fl  and  f’l  are  partial  derivatives  In  respect  to  y  and  y'  from 

V  Jr 

the  Integrand  f . 

If  the  Integrand  In  (7)  does  not  depend  on  x  *  the  Euler's 
equation  Is  Integrated  once  and  assumes  the  form 

fCy*  y')  *  f'fyify.  y*)  -  c,  • 

where  C  Is  the  Integration  constant. 

Applying  this  condition  to  Expression  (6)  we  easily  obtain 

[(R*  -  1)/rS1p(R')  -  C, 

where  the  designations  are  that  P(R')  *  (R')*9(R')  and  9'(R')  Is  the 
derivative  of  the  function  of  <p  given  by  the  graph  In  Pig.  1  In 
respect  to  R'  ■  tan  0/2. 

To  solve  Eq.  (8)  It  is  necessary  to  figure  out  P(R')  In  advance. 
Since  function  9(R')  Is  prescribed  by  the  graph*  the  values  of  derlva* 
tlve  7'(R')  needed  to  compute  P(R')  are  determined  numerically  with 
the  five  point  foxmila 

where  h  Is  the  range  of  values  of  R'  and  . . .  are 

respectively  the  first,  second,  etc.  differences.  The  function  P(R') 
obtained  from  the  computations  Is  shown  In  Fig. 

A  sufficient  minimum  condition  for  Integral  (7)  Is  the  fulflllMnt 
of  the  Inequality  ^  ^  along  the  extremal  y  *■  y(x)  . 

For  Expression  (6)  this  condition  has  the  forou 


1  1 

But  seeing  that  4  j^a  IT*  >  0>  fulfillment  of  the  last  condltlcm 
Is  equivalent  to  fulfilling  the  Inequality 


(9) 


Inequality  (9)  Is  fulfilled  only  on  the  ascending  side  of  cux*ve 

P(R')  (Plg.  3).  On  the  descending  side  of  curve  P(R')  we  have 
dP(R') 


dR' 


4!  0.  Consequently,  at  values  of  R'  >  where  R*^^^ 


Is  a  root  of  the  equation 

^  -0.  (10) 

no  such  smooth  curve  R(x)  exists  along  which  Integral  (6)  assumes  a 
minimum  value. 

It  Is  also  easy  to  see  that  Eq.  (8)  may  be  solved  beginning  only 
at  values  of  R  ^  Ro,  where  the  magnitude  Ro  depends  on  the  choice 
of  arbitrary  constant  C  entering  Eq.  (8) .  Actually,  by  rewriting 
Eq.  (8)  In  the  form 

It  Is  easy  to  satisfy  oneself  that  when  R-»  1  the  right  side  of 
the  resulting  equation  Increases  Infinitely  and  the  left  side  Is 
always  bounded  (Pig.  3).  Therefore  with  values  of  R  sufficiently 
near  to  1  Eq.  (8')  has  no  solution  and  consequently  with  the  Indicated 
values  of  R  there  Is  no  smooth  curve  realizing  the  minimum  of  Integral 
(6). 

Prom  what  has  been  set  forth  It  follows  that  In  the  class  of 
smooth  curves  of  R(x)  (to  be  differentiated)  there  Is  none  which  at 
the  prescribed  length  and  degree  of  expansion  gives  a  diffuser  wall 
profile  with  minimum  expansion  losses.  In  this  connection  It  Is 
natural  to  suppose  that  the  diffuser  wall  profile  with  minimum  losses 


does  not  represent  a  smooth  curve  but  forms  on  Intake  ( small  values 
of  R)  and  on  exit  ( large  values  of  B)  abrupt  area  changes  ( sudden 
expansion) ,  as  Is  shown  In  Fig.  4.  In  this  case  the  expression  for 
the  loss  factor  will  take  the  form: 


Here  the  first  and  last  members  represent  losses  In  sudden 
expansion  on  entry  and  exit,  respectively,  referred  to  the  kinetic 
energy  of  the  stream  on  entering  the  diffuser  (Fig.  4  shows  values 
of  Ro  and  rJ  .  If  the  curve  of  R  »  r(x)  when  0<  x<  Lis  so  chosen 

that  It  realizes  the  minimum  of  Integral 

t  I 

/-4 

then  coefficient  C  will  be  a  function  only  of  R©  and  R_  will  be  the 
ordinate  of  the  points  through  which  the  extremal  passes. 

The  values  of  Ro  and  R^  must  be  chosen  from  the  conditions  In 
Lavrent'ev  and  lyustemlk  (5lt 

^  -V  ^  « 

The  application  of  these  conditions  to  (ll)  gives: 


?'{«;>  4B+f(«i)- 1 -oj 


(12) 


Function  P(R')  on  the  left  side  of  Eq.  (8)  may  with  sufficient 
accuracy  be  approximated  by  a  power  relation  (Fig.  ?)  within  the 
limits  from  R*  -  0  to  a  value  of  R'^,^  determined  by  the  solution 
to  Eq.  ( 10)  . 


This  function  may  be  written,  for  example,  in  the  form 

F(R’)  -  51  (»’)••■• 


Then  Euler's  equation  for  the  Indicated  values  of  R*  will  have  t)w 


form: 


w(»y^-CjirrT-  (13) 

Hence  after  Integration  and  slight  transformations  we  obtain 


where 


<(R) 


dH 


Flguz^  5  shows  function  lO  . 

If  we  exclude  the  constant  C  from  Eq.  ( l4)  ( C  can  be  determined 
by  setting  x  -  L  and  R  »  R_  in  the  equation)  ,  we  will  have: 

O 


jt  ♦(«-♦(*•) 


(15) 


The  tangent  of  the  slope  angle  R*  nay  be  determined  from  Eq.  ( 1?) 
by  first  substituting  In  It  the  value  of  the  constant: 


\vm 


(16) 


- r 

Equations  ( 13)  and  ( l6)  determine  the  form  of  the  diffuser 
having  minimum  expansion  losses. 

In  order  that  these  equations  may  be  used,  the  method  of 
determining  the  magnitudes  of  Ro  and  R_  must  be  Indicated.  Let  us 
write  Eq.  (l6)  for  the  beginning  and  end  of  the  diffuser: 

•k-U-V  - ■ 

If  we  adduce  Eq.  ( 12)  «  Introduce  the  deslgnaticms 


(17) 


and  divide  the  first  of  Eqs.  ( 17)  by  the  aecond,  then  the  system 


( 12} ,  ( 17)  may  be  written  In  the  form 


(18) 


The  graphs  of  the  functions  4>(r)  and  A(  R')  are  shown  In 
Figures  6  and  3,  respectively.  Prom  the  first  equation  we  easily 
find  that  R'  =  0.177.  Further,  for  each  pair  of  values  n  and  R_ 

O  B 

from  the  second  equation  we  find  R^,  and  then  from  the  third  and 
fourth,  the  corresponding  values  of  Ro  and  L.  Prom  the  data  obtained 
we  construct  the  relations  Ro  =  fi(L)  and  R_  ■  fs(l)  for  several 
values  of  n  (Pigs.  7  and  8)  .  Let  us  remember  that  the  magnitudes 
Ro,  Rg,  and  L  refer  to  radius  R^  of  the  entry  section  of  the  diffuser. 

Prom  Figure  7  it  follows  that  the  value  of  Ro  for  the  usually 
used  diffusers  Is  practically  no  different  from  xailty.  Ihls  means 
that  the  sudden  expetnslon  on  entry  is  slight  and  can  be  neglected. 

As  for  the  sudden  expansion  on  exit,  that  Is  very  considerable  at 
small  lengths  (Plg.  8)  tends  to  zero  when  L  increases,  so  that  each 
curve  of  R^  =  L)  has  a  horizontal  asymptote  R^  ■  ^/TT.  Prom  Fig. 8 
It  also  follows  that  at  a  given  dimensionless  length  the  role  of 
sudden  expansion  Increases  with  the  Increase  in  n. 

Optimum  diffuser  shape,  as  was  shown  above,  is  In  the  part  where 
there  is  no  abrupt  expansion  determined  by  Eqs.  ( 15}  and  ( iQ . 
Notwithstanding  that  the  local  expansion  angle  determined  by  Eq.  ( 
is  a  variable  magnitude,  function  R(x)  Is  close- to  the  linear. 
Actually,  the  nature  of  function  rC^O  Is  determined  by  the  form 
of  function  if(B)  which  differs  very  sll|$itV  f rom  the  linear  (Fig.  5): 
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nius  the  optlaum  le  a  stepped  diffuser  consisting  of  a  part 
little  different  from  the  conical  with  a  degree  of  expansion  of 
R*  and  abrupt  expansion  on  exit. 

B 

The  work  of  Idel'chlk  [2]  shows  the  result  of  an  experimental 
ccniq>arlson  of  flat  diffusers  of  various  shapes*.  Of  the  three 
diffusers.  Isograde,  rectilinear  (6  •  const),  and  stepped,  with 
the  saaie  length  and  degree  of  expansion,  the  best  proved  to  be  the 
stepped  one.  The  same  work  displays  a  systematic  Investigation  of 
flat  Isograde  and  rectilinear  diffusers.  A  comparison  of  the 
efficiency  of  these  diffusers  shows  that  at  small  equivalent 
expansion  angles 

t--?— I- 

losses  In  the  rectilinear  diffuser  are  found  to  be  less  than  In  the 
Isograde.  With  substantial  equivalent  angles  the  Isograde  diffuser 
proves  to  be  BK>re  efficient.  These  conclusions  agree  with  the  results 
of  this  work.  Actually,  at  small  equivalent  angles  ( small  values 
of  n  and  large  ones  of  the  role  of  abrupt  expansion  at  the  exit 
of  an  optlaum  diffuser,  as  follows  fram  Fig.  8,  Is  small;  and  there¬ 
fore  the  rectilinear  diffuser  Is  nearer  the  optimum  than  the  Isograde 
Is.  Contrariwise,  with  large  values  of  0^^  (large  n  and  small 
abrupt  expansion  on  exit  from  the  optimum  diffuser  Is  found  to  be 
substantial  and  consequently  In  this  case  the  Isograde  diffuser  will 
be  nearer  the  optimum  than  the  rectilinear. 

In  conclusion,  let  us  note  that  In  the  examined  problem  only 
expansion  losses  were  taken  Into  consideration.  It  Is  easy, 

*Unfortunately,  similar  data  for  axlsynmetrlcal  diffusers  are 
lacking  in  the  literature  known  to  us. 

FT1>-TT-62-1738/1-^2«4 

_ _ _  _ _ HI 


however,  to  show  that  allowance  for  friction  losses,  which  In  the 
majority  of  cases  are  small  In  comparison  with  expansion  losses, 
does  not  practically  change  the  results. 

Conclusions 

The  variational  problem  of  the  shape  of  a  diffuser  with  minimum 
losses  has  been  examined  In  this  work. 

We  have  shown  that  the  optimum  Is  a  diffuser  consisting  of  a 
part  close  to  the  conical  and  an  abrupt  expansion  on  exit. 

The  results  we  obtained  agree  qualitatively  with  the  published 
experimental  data  on  flat  diffusers. 
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NOMOGRAMS  FOR  THE  APPROXIMATE  CALCULATION 
OF  PROFILE  FLOW  IN  TURBINE  GRILLS 

G.  Ya.  Shnltser 

A.  N.  Sherstyuk  [1]  has  proposed  an  approximate  method  of 
calculating  the  distribution  of  velocities  on  the  profile  perimeter 
of  turbine  grid  vanes  based  on  simple  solutions  for  the  flow  of 
an  incompressible  fluid  in  curvilinear  channels.  To  compute 
velocity  Wi  at  point  A  on  the  back  of  a  vane  (Fig.  l)  Sherstyuk  pro¬ 
poses  using  the  following  formula  after  having  inscribed  in  the  inter¬ 
vene  channel  a  circle  touching  the  vanes  at  A  and  B  and  after  having 
determined  the  width  of  the  channel  along  the  equlpotentlal  line 
AB: 

( 1) 

1  +  0.291  2(2  +  Sy-l^ 

Where  Wi  «  WiAf-i  ■  Ra/R»i  ^  ■  VRiJ  ^  ■  t/Rj;  Wj  is  the  stream 
z 

velocity  on  the  back  of  the  vane  at  point  A;  w,,  the  expenditure 

z 

component  of  velocity;  Ra«  the  radius  of  curvature  of  the  concave 
side  of  the  vane;  Ri,  the  radius  of  curvature  of  the  back  of  the 
vane;  h,  the  width  of  the  channel  measured  along  the  equlpotentlal  line 
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stream  velocity  w*  on  the  concave  side  of  the  vane  Is  computed 


fi*om  this  formulai 


1-0.8^ 

Wg  _ 

^  14-0.87 


(?) 


Despite  the  fact  that  two  works  [2«  3]  later  proposed 
somewhat  more  accurate  formulas  for  computing  flow  In  cui*vlllnear 
channels,  Sherstyuk's  method  Is  widely  used  also  at  the  present  tlsie 
for  designing  turbine  vanes  because  of  Its  simplicity  and  sufficient 
accuracy  for  qualitative  estimates.  It  Is  especially  useful  In 
cases  where  it  Is  requisite  quickly  to  calculate  the  natiire  of  the 
profile  to  be  used  and  also  for  obtaining  a  first  approximation  for 
more  rigorous  calculations. 

Therefore  It  Is  of  interest  for  shortening  the  computations 
to  determine  stream  velocities  fz*om  the  nomograms  proposed  below 
instead  of  f 3?om  Formulas  ( l)  and  ( 2)  . 

Let  us  rewrite  Formula  ( l)  In  the  form 


where 


I 


l  +  07Si  3(S-|-1^C 


(4) 


[nstead  of  ( 4)  we  may  write 


o.s«ll 


(5) 


The  last  expression  Is  an  equation  of  the  fourth  n<xiiographlo  order 
(in  the  Cauchy  fons)  in  respect  to  the  e3q>resslons  in  brackets. 

Such  equations  may  be  represented  In  the  form  of  nomograms  of  aligned 
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points 


On  the  basis  of  our  experience  with  the  calculations  we  will 
assume  the  following  limits  of  change  In  the  variables  In  Eq.  (3)t 

o.a<I<s.I 


In  addition^  the  following  may  be  encountered: 


Kt>moai  1 

o<*<oa/ 


(7) 


In  this  case  It  Is  more  convenient  to  use  Formulas  ( l)*  and  (2) 
which  are  considerably  simplified  and  take  on  the  form 


(8) 

(9) 


The  nomogram  of  Eq,  (5)  for  computing  the  value  of  Pj  with  re¬ 
spect  to  Ha  and  Ti  consists  of  the  rectilinear  scales  and  Ft 
Intersecting  each  other  at  a  right  angle  and  a  curvilinear  scale 
Hi  .  The  scale  markings  corresponding  to  the  values  of  He,  H  ,  and 
Pi  connected  by  Formula  (4)  lie  on  a  straight  line  (Plg.  2).  The 
scales  of  Hs  and  Ft  are  Isometric. 

The  curve  on  which  scale  hi  Is  located  Is  constructed  from 
equations 

,-T+o* 


(10) 


The  value  of  H  changes  from  0.2  to  2.0;  x  Is  plotted  along 
axis  Ft  and  £  along  axis  He  In  the  imlts  of  those  scales.  At 
points  X,  2.  ctunre  Hi  are  Inscribed  the  valtws  of  H  at  which  these 
points  were  cosiputed. 
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The  construction  of  the  Tii  scale  will  be  simplified  If  we  use 
the  hyperbolic  function  x  •  ( l/Ti)  -i*  0.29*  By  constructing  the  graph 
of  this  auxiliary  hyperbola  It  Is  easy  to  detenslne  the  value  of  H 
for  any  point  on  curve  Tii  as  the  ordinate  of  the  hyperbola  above 
this  point  (Fig.  2)  .  Thus,  on  the  Tit  scale  It  Is  sufficient  to 
write  down  only  basic  reference  points  (0.2,  0.?,  0.4,  etc.). 

After  the  value  of  Ft  has  been  found,  velocity  Wt  Is  determined 
from  Eq.  (3)  or,  more  conveniently  from  the  z-nomogram  coirrespondlng 
to  this  equation.  The  z-nomogram  consists  of  the  parallel  recti¬ 
linear  scales  Ft  and  Wt  and  the  rectilinear  scale  T  Intersecting 
them.  The  Ft  and  Wi  scales  coincide  with  the  Ft  and  Ra  scales 
of  the  preceding  nomogram.  The  t  scale  Is  projective  and  Is 
plotted  directly  from  the  eqtiatlon  “t  »  Wi/Ft.  The  scale  max4clngs 
corresponding  to  the  values  of  Ft,  T,  Wt  are  connected  by  Formula 
(3)  and  lie  on  a  straight  line.  For  convenience  In  use  both  nomograms 
are  combined  so  that  their  Ft  scales  coincide  and  may  be  blank. 

In  order  to  design  nomograms  for  calculating  the  velocity  on 
the  concave  side  of  the  vane  we  rewrite  Formula  (2)  In  the  form 


(11) 


where 


F, 


i+o.5r 


(12) 


Instead  of  ( 12)  we  may  write 

[7;li-"+[-2|rl['+ij+!-‘  (15) 

Equation  ( 13)  may  be  In  the  form  of  a  nomogram  of  aligned 
points  with  two  Intersecting  scales.  The  third  scale  in  this -case 
will  also  be  rectilinear  and  Isometric.  It  Is  constnieted  according 
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to  the  equations 


jr-.4  +  '*  (1*^) 

As  before  also,  x  Is  plotted  along  the  ?•  axis  and  £  along  the 
Ra  axis  In  units  of  these  scales. 

The  rule  for  using  the  nomogram  Is  evident  from  Fig.  2. 

After  the  value  of  Fa  has  been  found,  velocity  Wa  Is  determined 
by  Eq.  ( ll)  or,  more  conveniently,  from  the  z-nomogram.  In  the 
latter  case  the  already  mentioned  z-nomogram  Is  used  for  figuring 
Wi  with  only  this  difference,  that  Instead  of  Fi  the  value  of  Fa 
Is  plotted  on  that  scale,  and  the  value  Wa  Instead  of  T  .  Thus, 
all  the  nomograms  described  above  are  combined  (Fig.  2).  They  may 
easily  be  constructed  since  almost  all  the  scales  are  Iscsnetrlc. 

The  nomograms  guarantee  velocity  computations  to  an  accuracy  of 
three  significant  digits;  greater  accuracy  Is  superfluous,  as  the 
Initial  value  of  h  In  this  calculation  Is  taken  from  a  drawing  with 
the  same  accuracy. 

An  example  of  computation: 

From  the  drawing  are  taken  the  following  data:  *•  1.78, 

H  -  0.45,  t  -  0.83. 

Determining  velocity  Wi  .  Key:  If*  —  TIi-{Fi)— 

We  Join  point  1.78  on  scale  Ra  with  point  0.45  on  scale  Hi 
with  a  straight  line  which  we  extend  to  Intersect  with  the  P  scale. 
Between  the  point  thus  obtained  and  point  0.83  on  scale  T  we  draw 
a  stalght  line  at  whose  Intersection  with  the  Wt  scale  we  read  the 
answer-wi  *■  2.22. 

Determining  velocity  Wa.  Key:  T[a-ha--(Fa)-Wa  *^1. 

Point  1.78  on  vertical  scale  Ha  and  point  0.45  on  scale  Ha 
are  connected  by  a  straight  line  projected  to  Intersect  scale  P. 
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Through  the  point  thus  obtained  and  point  2.22  on  scale  Wi  we  draw 
a  straight  line  at  whose  Intersection  with  scale  Wa  we  read  the  an¬ 
swer  -  w«  -  1.58. 
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